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ABSTRACT 



We report on a timing analysis performed on a 62-ks long XMM-Newton observation of the accreting millisecond pulsar 
SAX J1808. 4-3658 during the latest X-ray outburst that started on September 21, 2008. By connecting the time of 
arrivals of the pulses observed during the XMM-Newton observation, we derived the best-fit orbital solution and a best- 
fit value of the spin period for the 2008 outburst. Comparing this new set of orbital parameters and, in particular, the 
value of the time of ascending-node passage with the orbital parameters derived for the previous four X-ray outbursts 
of SAX J1808. 4-3658 observed by the PCA onboard RXTE, we find an updated value of the orbital period derivative, 
which turns out to be P or b = (3.89 ± 0.15) x 10 _1 s/s. This new value of the orbital period derivative agrees with the 
previously reported value, demonstrating that the orbital period derivative in this source has remained stable over the 
past ten years. Although this timespan is not sufficient yet for confirming the secular evolution of the system, we again 
propose an explanation of this behavior in terms of a highly non- conservative mass transfer in this system, where the 
accreted mass (as derived from the X-ray luminosity during outbursts) accounts for a mere 1% of the mass lost by the 
companion. 

Key words, stars: neutron — stars: magnetic fields — pulsars: general — pulsars: individual: SAX J1808. 4-3658 - 
X-ray: binaries — X-ray: pulsars 



1. Introduction 

Accreting millisecond pulsars (AMSP) are generally inter- 
preted as the evolutionary link between low-mass X-ray 
binaries and rotation-powered millisecond radio pulsars. 
According to the recycling scenario, the last are in fact 
formed after a phase of mass accretion from a low-mass 
companion star, w hich ultimately sp ins the neutron star 
up to ms periods (jAlpar et al l 11982). SAX J1808.4-3658 
was t he first AMSP discovered (|Wiinands &: van der Klisl 
1998). In ten years the population of AMSP has grown to 
ten sources, all residing in close and transient binaries, but 
SAX J1808. 4-3658 can still be considered the cornerstone 
of its class as it has repeatedly gone in outburst almost ev- 
ery two years after its first detection, making it the most 
obscrvationally rich source of its class. Its timing and or- 
bital properties have been extensively studied; while the 
timing behavior of the spin frequency is still debated be- 
cause of the pr esence of timing noise i n the phase delays (see 
IBurderi et atll2006l B06 hereinafter: iHartman et~aUl2008l 
H08), there is good agreement as regards the orbital pa- 
rameters of this source, which are no w known with extreme 
precision (H08: lDi Salvo et al.li2008l . DS08). 

In particular, a timing of the past four outbursts of SAX 
J1808.4-3658 observed with the PCA onboard the RXTE 
satellite (spanning more than 7 years from April 1998 to 
October 2005) has allowed the orbital period of this source 
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to be derived with relative uncertainty AP or (,/P or (, of ~ 
2 x 10" 9 (dS08, H08). More interesting, the orbital period 
of the source shows a significant derivative, which indicates 
that it is increasing at a large, as well as puzzling, rate of 
-Porb = (3.5 ± 0.2) x 10~ 12 s/s. The value of the orbital 
period derivative is at least one order of magnitude greater 
than what is expected for conservative mass transfer driven 
by gravitational radiation (GR; DS08, H08). 

To continue the monitoring of the orbital period deriva- 
tive and to further improve the orbital solution for this in- 
teresting system, we analyzed a 62 ks XMM-Newton Target 
of Opportunity (ToO) observation of this source, performed 
during the latest outburst in October 2008. 



2. Timing analysis and results 

SAX J1808.4-3658 was found in outburst on 2008 
September 21 by RXTE , and since then it has been the ob- 
ject of an intensive observational campaign. XMM-Newton 
observed SAX J1808.4-3658 as a ToO observation for 62 ks 
on 2008 October 1 (start time MJD 54739.99517), roughly 
one week after the assumed outburst peak. In this paper we 
concentrate on the pn data, which have the best statistics 
for timing analysis. A detailed description of the XMM- 
Newton obs ervation and data r eduction has already been 
reported in iPapitto et al.l (|2009D . The arrival times of all 
the pn events were referred to the solar system barycenter 
by using the barycen tool in SAS v. 8. 0.0. As the best esti- 
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mate for the source coordinates, we considered the position 
of the radio co unterpart, which ha s an uncertainty of 0.4 
arcsec (90% c.L. lRupen et a l. 2002) and is compatible with 
the optical counterpart ( Giles et al.lll999n . 

We used the XMM-Newton/pn data to improve the or- 
bital solution for SAX J1808. 4-3658. To this aim we also re- 
analyzed the previous four outbursts of this source observed 
by RXTE/PCA, namely 1998 April (Obs. ID P30411), 
2000 February (Obs. ID P40035), 2002 October (Obs. ID 
P70080), and 2005 June (Obs. ID P91056 and Obs. ID 
P91418) outbursts, respectively (see H08 for a detailed de- 
scription of these observations). With this new analysis we 
present a more conservative treatment of the errors on the 
phase delays, to avoid any possible underestimation of the 
uncertainties on the orbital parameters caused by the phase 
changes and shifts that are known to occur during out- 
bursts. 

For each outburst, we repeat the analysis performed 
by DS08, with a few differences. Differential corrections to 
the adopted orbital parameters were found by fitting the 
pulse phases as a function of time for each outburst (see 
e.g. Fig[T]). In general, any residual orbital modulation is 
superposed onto a long-term variation in the phases, e.g. 
caused by a variation of the spin. However, SAX J1808.4- 
3658 shows very complex behavior of the pulse phases with 
time, with phase shifts, probably caused by variations in 
the pulse shape, that are difficult to interpret (see e.g. B06, 
H08). To model this complex behavior of the pulse phases 
with time, we fit the phases of each of these intervals with 
the formula for th e differential correc tions to the orbital pa- 
rameters (see e.g. lDeeter et al.l ll981,) plus a polynomial (up 
to 8 th de gree, depending on the irregular behavior of the 
pulse phases) . This technique is complementary to the tech- 
nique adopted in dS08, where we preferred to restrict the 
fit of the differential corrections to the orbital parameters 
to intervals in which the long-term variation and/or phase 
shifts were negligible and considered consecutive intervals 
with a duration of at least 4 orbital periods (depending 
on the statistics), fitting the phases of each of these inter- 
vals only with the formula for the differential corrections 
of the orbital parameters. We have verified that the two 
techniques give similar orbital corrections, but we now ob- 
tain somewhat larger uncertainties on the derived orbital 
parameters. 

Over the 10-year base line SAX J1808.4-3658 has been 
observed, the measured P or b should have produced ~ 1 ms 
change in P or b- However, the uncertainty on the individual 
P or b determinations from the local orbital solution during 
each outburst are about 0.5-2.0 ms (see also H08); accord- 
ingly, no significant corrections to P or b have been directly 
observed. Thus we fix our parameter P or b,o = 7249.156499 
s (dS08) and improve its measure below using the same 
technique as in dS08. For each of the five outbursts, we 
find similar corrections for the projected semimajor axis of 
the neutron star orbit and upper limits on the eccentricity. 
We hence combine the measurements corresponding to each 
of the outbursts by computing the error-weighted means of 
the corrections to the projected semimajor axis and of the 
upper limits on the eccentricity. On the other hand, we find 
that the correction to the predicted times of passage of the 
NS at the ascending node at the beginning of each outburst, 
T m predicted = To+NP olh Q , is different among the five out- 
bursts. Here is the time of ascending node passage given 
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Fig. 1. Phase delays (in fis) vs. time during the XMM- 
Newton observation after the first iteration. To show for 
clarity the orbital correction against an horizontal line, the 
frequency adopted for the folding is the best-fit spin fre- 
quency of the 2008 outburst. The orbital modulation is 
clearly visible (top panel) and is well-fitted by our orbital 
solution, as can be seen in the post-fit residuals (middle 
and bottom panels). 



by dS08 Q, and the integer N is the exact number of orbital 
cycles elapsed between two different ascending node pas- 
sages; i.e., N is the closest integer to (T^— P *)/P or b o under 
the assumption that \T* n - T* n prcdictcd | << Porb o that we 
have verified a posteriori, and m = 1998,2000,2002,2005, 
and 2008. These points show a clear parabolic trend that 
we fitted to the formula 



8T r * n = 5T* + SP oth xN+ (l/2)P orb Porb o x A^ 2 



(1) 



where 5T*, (5P or b, and P or b (the orbital period derivative) 
are the fit parameters. To compute the correct errors on the 
derived parameters, we also considered the errors induced 
by the uncertainties in the time coordinate (which we as- 
sumed to be half of the duration of each outburst) following 
the stan dard procedure for the propa gation of these errors 
(see e.g. iBevington fc Robinsonll2003h . 

We iterate the method described above until the upper 
limit on the eccentricity, as well as the corrections on the 
projected semimajor axis, the orbital period, and its deriva- 
tive were all compatible with zero within the computed er- 
rors. In Fig. Q] we plot the phase delays obtained for the 
XMM-Newton/pn data with the orbital solution of DS08 
neglecting the orbital-period derivative. The sinusoidal os- 
cillation is clearly visible (top panel) and mainly stems from 
the ~ 35 s delay in the ascending node passage with respect 
to the one predicted with a constant orbital period. In the 
middle panel, wc show the time residuals after our first im- 
provement of the orbital solution (reduced x 2 — 1-86 for 49 
d.o.f.). In the bottom panel, the time residuals are shown 
with respect to the final orbital solution shown in Table [T] 
(reduced \ 2 — 1-07 f° r ?2 d.o.f.). In this case more points 
are visible in the residuals. This is caused by how improv- 
ing the orbital solution results in a higher statistical signif- 



1 The actual T * used in this work is the T * of dS08 decre- 
mented by Porb o in order to have the time of ascending node 
passage just before the beginning of RXTE data of the 1998 
outburst which occurred at To = 50914.8099 MJD. 
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Table 1. Best-fit orbital solution for SAX J1808.4-3658 
derived from the analysis of the five outbursts observed 
from 1998 to 2008. 
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Fig. 2. Differential correction to the time of passage of the 
neutron star at the ascending node for each of the five out- 
bursts analyzed (top panel) and residuals with respect to 
the best fit parabola (bottom panel). 



icance for each point, making acceptable phase-points that 
were rejected with our previous solution. In Fig. [5] we plot 
the final ST^s vs time, together with the best-fit parabola. 
The fit was good with a reduced \ 2 = 1-55 (f° r 2 d.o.f.). 
This corresponds to a probability of 21% of obtaining a \ 2 
that is larger than the one we found. Our result is therefore 
acceptable, since the probability we obtained is well above 
the conventionally accepted s ignificance level of 5% (see e.g. 
iBevington fe Robinsonll2003h . 

In agreement with previous results, we find a highly 
significant derivative of the orbital period, which indicates 
that the orbital period in this system is increasing at a rate 
ofP orb = (3.89±0.15)x 10" 12 s/s. The best-fit values for the 
orbital parameters are shown in Table [TJ With these new 
ephemerides, we correct the events of the XMM-Newton 
ToO observation of 2008 October 1 and fit the 700 s folded 
phase delays of the fundamental with a straight line to ob- 
tain the correction to the spin frequency adopted for the 
folding. The fit is acceptable with a reduced x 2 = 1.16 (for 
71 d.o.f.). The best-fit spin frequency of the 2008 outburst 
is also reported in Table [TJ 



3. Discussion 

In this paper we present a timing analysis of the XMM- 
Newton observation performed during the last outburst 
of the AMSP SAX J1808.4-3658, together with a re- 
analysis of the past four outbursts from this source ob- 
served by RXTE/PCA. We find that the measures of the 
time of ascending node passage for the last 5 outbursts 
from this source are again perfectly fitted by a parabola; 
from this fitting, we find a derivative of the orbital period, 
-Porb = (3.89 ± 0.15) x 10~ 12 s/s, w hich agrees wi t h pre- 
viou s measures (H0 8 , DS0 8; see also iPatruno et all 



and lHartman et al.l (|2009f ) , who report a similar value us- 
ing RXTE data taken during the 2008 outburst). The con- 
clusion from these measurements is that the orbital-period 
derivative in this system has been stable over the past 
ten years. This supports the hypothesis that the system is 
following a non-conservative secular binary evolution (see 
DS08). 



Parameter Units 



Value 



2?998 


MJD 


50914.79452952(85) 


Porb 


s 


7249.156444(23) 


Porb 


1CT 12 s/s 


3.89(15) 


ai sin i/c 


s 


0.0628106(20) 
< 1.2 x 10~ 4 


e 




^2008 


Hz 


400.9752084532(17) 


P2OO8 


MJD 


54739.98342965 



Note: Errors are at 90% c.l. on the last 2 digits. The upper limit 
on the eccentricity is at 95% c.l. The values of P or b and P or b are 
referred to Ti 99S , while the value of ^2008 is referred to T2008 ' 



H08 note that the orbital parameters of SAX J1808.4- 
3658 are very similar to those of the so-called black-widow 
milli second pulsars, which show large and variable P or b (see 
e.g. lArzoumanian et"aT1 11994 iDoroshenko et all 1200 1[ ). In 
those cases P or b changes have been ascribed to gravita- 
tional quadrupole coupling, i.e. a variable quadrupole mo- 
ment of the companion that is caused by a cyclic spin- 
up and spin-down in the u pper layers of the companion 
(|Applegate fe Shahamlll99l . In this scenario, the compan- 
ion star must be partially non-degenerate, convective, and 
magnetically active, so that the wind of the companion star 
will result in a strong torque, which tends to slow the star 
down. This torque is then transferred to the orbit via grav- 
itational interaction. Typical timescales for the orbital pe- 
riod (roughly sin usoidal) modulation are about 5 — 6 yr 
(|Applegatei[l992T ). However, SAX J1808.4-3658 has been 
monitored since its discovery as a millisecond pulsar in 1998 
and does not show any evidence of this complex behavior 
in its orbital period. Instead, its P or b has been remarkably 
stable over this time. We cannot exclude, however, that its 
Porb might exhibit more complex behavior on decade-long 
time scales. 

The value of the orbital period derivative of SAX 
J1808. 4-3658 clearly indicates that the binary system is 
expanding at a rate that is more than one order of mag- 
nitude higher than what is predicted by conservative mass 
transfer driven by GR (H08, DS08). In fact, the angular 
momentum losses caused by GR have the effect of reduc- 
ing the binary orbital period. The transfer of mass from 
the secondary to the neutron star indeed has the effect of 
expanding the system (si nce the second ary star is lighter 
than the neutron star, e.g. IVerbuntJll993D . Following DS08, 
angular momentum conservation together with the third 
Kepler's law and considering the GR angular momentum 
losses, give the orbital period derivative: 



P_ 12 = 0.138 mfVi(l + 0.1q .ir 1/3 P 2 h 5/3 + 
6.840 m^ 1 qQ 1 1 P 2 h g(/3,q,a) to_ 9 



(2) 



where P_i2 is P or b in units of 10 12 s/s, mi = Mi/M Q , qo.i 
is the mass ratio q = TO2/TO1 in units of 0.1, P21, = P or b/2h, 
m_g = — M2/(lO -9 M0/yr), (3 is the fraction of the mass 
lost by the companion that is accreting onto the neutron 
star, a the specific angular momentum of the mass lost by 
the system written in units of the specific angular momen- 
tum of the companion (a = [1 - 0.462(1 + q) 2 / 3 ^ 1 / 3 ] 2 ~ 0.7 
for matter leaving the system with the specific angular mo- 
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Fig. 3. Mass-loss rate of the companion of SAX J1808.4- 
3658 in conservative (J3 = 1) and highly non-conservative 
(f3 = 0.01) scenarios. The gray area on the left is forbid- 
den because the radius of the (Brown Dwarf) companion 
is greater than its Roche Lobe. The gray area on the right 
shows the constraints on the companion mass imposed by 
the mass function of the system. 



mentum of the inner Lagrangian point), and g((3,q,a) = 
1 — 0q — (1 — 13) {a + q/3)/(l + q). Inserting our derived 
value P-12 = 3.89 into this equation, we can derive m_g vs. 
rri2 — qvn,\. The corresponding curves (very insensitive to 
the adopted values m\ — 1.1, 1.56, 2.2) are shown in Fig. [3] 
and labeled with [3—1 for the conservative mass transfer 
case. In the same figure the horizontal line at m_g = 0.017 
is the time-averaged mass accretion rate derived from the 
X-ray luminosity observed during outbursts and quiescence 
periods in the last ten years (which could even be an over- 
estimate by a factor 2 of the actual mass accretion rate, cf. 
iGallowav & Cummind l200l . The crossing point of these 
curves, which is the necessary solution imposed by the con- 
servation of angular momentum, falls in the forbidden re- 
gion at m 2 ~ 0.003. 

We can further impose the obvious constraint R2 < 
Rrl 2 > where R2 is the radius of the companion and i?RL 2 = 

1.2 x 10 10 m^ /3 P 2 2 h /3 cm is the radius of the Roche Lobe of 
the companion. Adopting results from detailed numerical 
calculations of the mass radius relation for low main se- 
quence stars and brown dw arfs of minimum radius (na mely 
cold, or older than 5 Gyr, IChabrier fc B araffc 2000), the 
gray area on the left in Fig. [^indicates the region excluded 
by the violation of this constraint and shows that a con- 
servative solution is impossible. Companion masses below 
0.02 M & are also excluded by the constraints imposed by 
the mass function of the system for reasonable values of the 
neutron star mass. 

In line with DS08, we thus propose a highly non- 
conservative scenario in which the mass-loss rate from the 
companion is stable at a level of m_g ~ 1, with accretion 
episodes lasting a few tens of days, separated by quiescent 
phases in which the same rate of mass is ejected by the sys- 
tem. In this hypothesis, we calculate that (3 — 9.5 x 10~ 3 . 
The corresponding lines (which are again computed for 
mi = 1.1, 1.56, 2.2) are shown in Fig. [3] and labeled with 
(3 = 0.01. The crossing point between the non-conservative 



lines and the horizontal line for m_g = 1.8 (from B06) at 
TO2 ~ 0.08 is outside the forbidden area and demonstrates 
that our proposed scenario is, in principle, viable. The ejec- 
tion of matter out of the system via a pulsar wind will not 
imply an extra spin down for the pulsar, since in our model 
matter is ejected by the radiation pressure of the pulsar, 
well outside the light cylinder radius (where the electro- 
magnetic radiation from the pulsar is no longer connected 
to the neutron star magnetic field). Also, any possible spin 
down caused by the propeller effect at the end of an out- 
burst may be averaged to zero by a spin up, because of the 
accretion inside the corotation radius, that is expected at 
the beginning of an outburst. 

In conclusion, to explain the large and stable orbital 
period derivative in the case of SAX J1808. 4-3658, we 
propose that mass is expelled at a high rate from the 
system with the specific angular momentum of the inner 
Lag rangian point (or, ev entually, of an outer disk, if present, 
e.g. iDelove et al.l [20081 . These results strengthen the hy- 
pothesis that SAX J1808. 4-3658 belongs to the population 
of the so-called 'hidden' millisecond pulsars, whose radio 
emission is completely blocked by material engulfing the 
system that is continuously replenished by the mass out- 
flow driven by companion irradiation. In the case of SAX 
J1808. 4-3658, irradiation is probably caused by the power 
emitted by the magncto-dipolc rotator, which may also ex- 
plain why, although the companion star is transferring mass 
at a high rate, this mass does not accrete onto the neutron 
star (see DS08 for a more detailed discussion) . If this is the 
case, ~ 99% of the transferred mass in this system is not 
directly observable. Although 10 years is not long enough 
to confirm that the system is in a stationary (secular) or- 
bital evolution, the next few years will be important for 
confirming the stability of the orbital period derivative or 
to observe P reversal as in classical black-widow pulsars. 

Acknowledgements. We thank the XMM-Newton team who sup- 
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